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ABSTRACT 

We use state-of-the-art numerical simulations to explore the observability and the 
expected physical properties of the progenitors of the local group galaxies at z ~ 6 — 8, 
within 1 billion years of the big bang. We find that the most massive progenitors of the 
Milky Way (MW) and Andromeda (M31) at z ~ 6 and 7 are predicted to have absolute 
UV continuum magnitudes Muv — —17 to —18, suggesting that their analogues lie 
close to the detection limits of the deepest near-infrared surveys conducted to date (i.e. 
HST WFC3/IR UDF12). This in turn confirms that the majority of currently known 
z ~ 6 — 8 galaxies are expected to be the seeds of present-day galaxies which are more 
massive than L* spirals. We also discuss the properties of the local-group progenitors 
at these early epochs, extending down to absolute magnitudes Mjjv — —13. The most 
massive MW/M31 progenitors at z ~ 7 have stellar masses ~ 10 8 Mq, stellar 
metallicities Z* ~ 3 — 6%Z Q , and predicted observed UV continuum slopes f3 ~ —2.4 
to -2.5. 
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1 INTRODUCTION 

In the last three years, the sensitive near-infrared imag- 
ing made possible by Wide Field Camera 3 (WFC3/IR) 
on the Hubble Space Telescope (HST) has tr ansformed the 
study of galaxies at z ~ 7 (see lDunloDll2012l ). Several hun- 
dred Lyman-break galaxies (LBGs) have now been uncov- 
ered with HST in the redshift range 6.5 < z < 8.5 (e.g . 



Oesch et al. 2010; Bouwens ct al. 2010a; McLure et al.120101: 



Finkelstein et al.ll201Cl; iMcLure et al.ll201ll: lLorenzoni et al 



20111 : iBouwens et al.ll201ll : IOesch et al.ll2012l : IMcLure et al 
2012J), and the new ground-b ased wide-field near-inf rared 
surveys, such as Ultra VISTA (|McCracken et all I2012T ) are 
now reaching the depths required to reveal galaxies more lu- 
minous than Muv — —19 at z — 7 (|Castellano et al.l l2010l : 
lOuchi et al.ll20ld : iBowler et al1l2012l ) . 

Most recently, attention has been refocused on pushing 
to even fainter magnitudes and still higher redshifts through 
the ultra-deep WFC3/IR imaging in the HUDF (GO 12498; 
hereafter UDF12). The recently-completed 128-orbit UDF12 



observations reach 5a detection limits of Yiog = 30.0, 
J125 = 29.5, J140 = 29.5, -Hi60 = 29.5 (after combination 
with the UDF09 da ta), and are the deepest near-infrared 
images ever t aken jEllis et all 120131: iDunlop et al.l l2ul2bl ; 
IMcLure et all 120121 : fechenker et al.ll2012l : lOno et al.lBoLW . 
A detailed description o f the UDF12 data-set is provided by 
iKoekemoer et al) (|2012l ). and the final reduced images are 
available on the UDF12 team web-page Q. 

The ever-improving data on galaxies at early cosmic 
times has motivated attempts to link these objects to the 
galaxies seen in the local Universe, and in particular to 
the Milky Way (MW) or indeed to L* galaxies in gen- 
eral (which dominate the local luminosity density). The 
approaches used to establish such a link range from the 
use of correlation functions to state-of-the-art cosmologi- 
cal simulations. For exa mple, traci n g me rger trees using 



cosmological simulations, iNagaming (|2002l ) found that the 
brightest LBGs (with My ^ —23) at z = 3 merge into 
groups and clusters at z = 0, and about half the LBGs 
with —23 ^ My ^ —22 evolve into local L* galaxies. Mean- 
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while, using a clusterin g analysis of a sam p le of 2 8500 galax- 
ies at z ~ 1.4 — 3.5, lAdelberger et all (|2005t ) concluded 
that the typical objects in their sample would have evolved 
into elliptical galaxies at z = 0. More recently, coupling 
N-body Dark Matter (DM) simulations with a prescrip- 
tion for star formation, feedba ck and chemical evolution, 
lOkrochkov fc Tumlinsonl (|2010l ) have claimed that the pro- 
genitors of the MW could be visible as LBGs up to z = 7 and 
2 = 9, including and excluding the effects of dust absorp- 
tion on ultra violet (UV) photons, r es pectively. The reade r 
is referred to ISalvadori et all d2010T). iGuaita et all J2010T) . 



Daval fc Libeskindl J2012T ). IWalker-Soler et al.l l|2012f ) and 



Yaiima et al.l (|2012) for calculations aimed at establishing 



whether local group (LG) progenitors could be detected as 
Lyman Alpha Emitters (LAEs) in existing or future surveys. 

In this study, our aim is to investigate the link between 
the faintest high-redshift LBGs uncovered to date (i.e. galax- 
ies with Muv ~ -17AB mag detected through the UDF12 
campaign) and local L* galaxies. Owing to computational 
limitations, at present, it is too costly to run a simulation 
with a large enough box size so as to be able to treat long- 
wavelength modes effectively, as well as a high enough spa- 
tial resolution to resolve galaxies as faint as those detected 
in the UDF12 imaging. We therefore employ the approach 
of using re-simulations of the LG: by construction, these re- 
produce the properties of the two local ~ L* galaxies, the 
MW and M31 (Andromeda) at z — 0, while the small vol- 
ume simulated allows a high resolution to be achieved, as is 
discussed in more detail in what follows. 



2 THE SIMULATIONS 

For the calculations presented in this work, we have used the 
CLUES si mulation, which was run with the PMTree-SPH code 
GADGET2 (|Springe]||2005r ) in a cosmological box of size 64/i _1 
comoving Mpc (cMpc). The runs used sta ndard ACDM 
initia l conditions and WMAP3 parameters (|Spergel et al.l 
l2007h such that Q m = 0.24, Q b = 0.042, Q A = 0.76, 
h = 0.73, cr 8 = 0.73 and n = 0.95. 

In brief, observations of objects in the local universe 
were used to reconstruct t he initial density field u sing the 
Hoffman- Ribak algorithm IjHoffman fe Ribald 119911 ). These 
constrained initial conditions force the 2 = linear scales 
to resemble the input, the local universe; non-linear scales 
are unconstrained. To obtain a reliable LG, a number of 
low-resolution runs were performed until a suitable LG (de- 
fined in terms of the mass, relative distance and velocity of 
the two MW and M31 mass haloes) was found. These ini- 
tial conditions were then re-run with high r esolution and 
as dy namics using the prescriptions given bv lKlypin et al.l 
2001), to produce two objects resembling the MW and 
M31, each of which contains about 10 6 parti cles within their 
virial ra dii; interested r eader s are referred to Libeskind et al.l 
j|2010h, iKnebe et al.l (|201Ch . iLibeskind et all (|201lh and 
iKnebe et all ( 2011 ) for complete details. Using a spatial 
resolution of 150 pc, the mass resolution obtained was 
2.1 x 10 5 /i _1 M Q and 4.4 x 10 4 /i -1 M Q for DM and gas par- 
ticles, respectively. 



The simulation emp loys the feedback rules of 
ISpringel fc Hernquistl |2003l ): hot ambient gas and cold gas 
clouds in pressure equilibrium form the two components 
of the interstellar medium (ISM). Gas properties are cal- 
culated assuming a uniform but evolving UV backgroun d 
generated by QSOs and AGNs (jHaardt fc Madad Il996r i. 
Metal line (or molecular) dependent cooling below 10 4 K is 
ignored. Star formation is treated stochastically, choosing 
model para meters that reproduce th e Kennicutt law for spi- 
ral galaxies (|Kennicuttlll983l .l 1998). The instantaneous re- 
cycling approximation is assumed: cold gas cloud formation 
(by thermal instability), star formation, the evaporation of 
gas clouds, and the heating of ambient gas by supernova 
driven winds all occur at the same instant. We assume ki- 
netic feedback in the form of winds driven by stellar explo- 
sions. 

Haloes and subhaloes are identifi ed using the publicly 
availa ble Amiga Halo Finder (AHF; iKnollmann fc Knebd 
2009). AHF uses an adaptive grid to find local over-densities; 
after calculating the gravitational potential, particles that 
are bound to the potential are assigned to the halo. 

In order to find the progenitors of the MW and M31 
that could be visible as LBGs at z ~ 6 — 8, we identify all 
the particles within their virial radii at z = 0. These parti- 
cles are then followed back in time and located at z ~ 6, 7 
and 8 in the simulation box; if bound to a structure, the cor- 
responding halo is considered to be the progenitor of a local 
L* galaxy at that redshift. Further, in all the calculations 
presented in this paper, we only use those high-redshift pro- 
genitors of the MW and M31 that contain at least 20 star 
particles; this corresponds to galaxies that contain a mini- 
mum of 540 total (DM, gas and star) particles at z ~ 6 — 8. 

Although we have used an extremely high-resolution 
simulation that reproduces the masses, relative velocities 
and separation of the two local L* galaxies (the MW and 
M31) to trace their progenitors back in time, the main limi- 
tation of this work is that all the results shown in this paper 
relate to the specific merger tree obtained from the simu- 
lation. It would be interesting to explore how these results 
might change in the future when/if a statistically-significant 
number of similar resolution LG simulations become avail- 
able. 



2.1 Identifying LBGs 

To identify the progenitors of the MW and M31 that could 
be observed as high-redshift LBGs, we start by computing 
their UV luminosities (or absolute magnitudes): we consider 
each star particle to form in a burst, after which it evolves 
passively. The total spectral energy distribution (SED), in- 
cluding both the stellar and nebular continuum, for each 
star particle i s computed via the p opulation synthesis code 
STARBURST99 (|Leitherer et alj|l999f ). using its mass, stellar 
metallicity and age. The age is computed as the time differ- 
ence between the redshift of the snapshot and the redshift of 
formation of the star particle. Further, we assume a zero es- 
cape fraction of neutral hydrogen ionizing photons, thereby 
maximizing the contribution of the nebular continuum to 
the SED. The total intrinsic continuum luminosity, L % " (at 
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Figure 1. The mass- weighted stellar metallicity, Z», as a function of the stellar mass, M„, for the simulated progenitor galaxies at 
z ~ 6 — 8. Open circles and triangles show the Z» — M* relation for the progenitors of the MW and M31 respectively. The points have 
been colour-coded according to the absolute UV magnitude shown on the colour-bar. 



1500 A in the galaxy rest frame) for each progenitor is then 
calculated by summing the SEDs of all its star particles. 

Continuum photons can be absorbed by dust within the 
ISM and only a fraction, f c , escape out of any progenitor, 
unattenuated by dust. The value of f c for each progenitor 
is calculated assuming Type II supernovae (SNII) to be the 
primary dust f actories such that each SNII produces about 
0.4M Q of dust ijTodini fc Ferrarall200lh . each SNII destroys 
dust with an efficiency of about 20% within the re gion it 
shocks to speeds > 100 kms" 1 (|Seab fc ShuJll983h . a ho- 
mogeneous mixture of gas and dust is astrated into star 
formation, and that dust is lost in SNII-powered outflows. 
We assume the dust to be distributed in a screen of ra- 
dius Vd = r g where the gas distribution radius is calculated 
to be r g = 4.5Ar2oo- Here, th e spin parameter has a value 
A = 0.05 IjFerrara et al.ll2000f ) . and r2oo is the virial radius. 
This dust distribution radius is used to obtain the dust sur- 
face mass density which can then be translated in the dust 
optical depth to obtain f c (for details of this calculation see 
iDaval et al.ll20ld . l2012h . The observed continuum luminos- 
ity L° bs can then be expressed as L° bs = L l c nt x f c . 



3 RESULTS 

As a result of the exquisite simulation resolution we are 
able to study LBGs that are nearly two orders of magni- 
tude fainter in luminosity (Muv < —13) than the current 
near-infrared detection limit of Muv — —17 at z ~ 7; such 
faint LBGs could potentially be detected with future facili- 
ties such as the James Webb Space Telescope (JWST). We 
now present the theoretical results for the physical prop- 
erties of the MW and M31 progenitors, and the resulting 
predictions for observables such as the luminosity functions 
(LFs) and UV spectral slopes, /3 (fx oc A* 3 ), for compari- 
son with the most re cent observational r esults presented by 
iMcLure et all (|2012h and lDunlop et al.1 l|2012bl ). 



3.1 The mass-metallicity-magnitude relation 

We start by showing the relation between the stellar mass 
(M»), mass- weighted stellar metallicity (Z*) and the dust- 
attenuated absolute UV magnitude (Muv) for the progen- 
itors of the MW and M31, in Fig. Q] As expected from the 
hierarchical structure formation scenario where successively 
larger systems assemble from the merger of smaller ones with 
time, the M* range increases with decreasing redshift, going 
from 10 5 - 7 " 7 5 M Q to lO 5 7 " 8 M Q in the 300 Myr between 
z ~ 8 and z ~ 6. This is accompanied by a significant in- 
crease in the number of identifiable MW and M31 progeni- 
tors, from about 64 at z ~ 8 to 210 at z ~ 6, as more and 
more tiny building blocks assemble with time and undergo 
star formation. Since metals are produced by star formation, 
progenitors that have built up a larger stellar mass are also 
more metal rich, giving rise to a mass-metallicity relation 
such that Z, increases from 10 -19 Z for M» = 1O 5 7 M 
to about 10~ 108 Z Q for M* = 1O 8 '°M at z ~ 6; the av- 
erage stellar metallicity rises from about 0.02 Zq at z ~ 8 
to 0.035 Z at z ~ 6. As a result of their larger potential 
wells (that prevent gas loss due to negative feedback from 
SN-driven winds), more massive galaxies contain larger gas 
masses available for star formation, and hence are the bright- 
est in the rest-frame UV. 

To summarize, we find a three-dimensional relation be- 
tween M* — Z, — Muv such that the more massive pro- 
genitors are both more metal enriched and brighter in the 
UV; these trends hold true for the progenitors of both the 
MW and M31. However, since the simulated stellar mass 
of M31 (1.98 x 1O 1O M ) is larger than that of the MW 
(1.62 x 1O 1O M ) at z = 0, the most massive progenitors 
belong to M31 at any of the redshifts studied. 

Because a few of the objects shown in Fig. 1 have 
Muv —17, we would predict that analogues of some of 
the most massive progenitors of the MW and M31 at z ~ 6, 
and of M31 at z ~ 7 are bright enough to be found among 
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Figure 2. The UV LFs of the progenitors of the MW and M31 at z ~ 6 — 8. In each panel, points show the observational results, while 
the dashed (solid) histograms show the intrinsic (dust atte nuated) UV LF s obta ined using the progenitors of both the MW and M31. 
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limits for HST (JWST) near-infrared imaging. 
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the faintest high-redshift LBGs detected by the HST at these 
redshifts. However, at z ~ 8, none of the progenitors of the 
two L* galaxies at z = are brighter than the current de- 
tection limit of Muv < —17.25, as is also shown in Fig. [2] 
that follows. 



3.2 UV luminosity functions 

In Fig.[2]we compare the predicted UV LF of the progenitors 
of the MW and M31 to the observed UV LF at z = 6 - 8. 
Within the re-simulated region, we find the centre of mass 
of the MW and M31 and build spheres that enclose 95% 
of the total mass assembled at that redshift: the combined 
volume occupied by the progenitors of the MW and M31 is 



then found to be « (120, 200, 175) cMpc 3 at 



(6,7,8). 



However, to estimate the contribution of the LF of the pro- 
genitors to the overall LF in a representative volume, we 
need to account for the typical comoving volume occupied 
by galaxies such as the MW and M31. 

At z = 0, the MW and M31 have ab solute if-band 
magn itudes of Mk = —24.02 and —24.51 l|Hammer et al.l 
|2007| ). To estimate the present-d ay number densi t y of s uch 
galaxies we use the jf-band LF of lKochanek et aL I (20011) af- 
ter correcting their results to the Hubble constant adopted 
in our simulation; a fter boosting the absolu te magnitudes 
given in Table 2 of iKochanek et al.l |200ll ) by 0.68 mag, 
and reducing the values of logn by —0.41 (i.e. to correct 
to h — 0.73), we find that the number density of galaxies 
in the absolute magnitude range —24.7 < Mk < —23.7 is 
~ 6.2 x 10 _3 cMpc -3 mag -1 . This implies that oneL* galaxy 
at z — typically occupies a volume of about 160 cMpc 3 , im- 
plying that the total volume occupied by the MW and M31 
is approximately 320 cMpc 3 . Normalizing the combined vol- 
ume occupied by the progenitors of both the MW and M31 
at z ~ 6, 7 and 8 to the volume of 320cMpc 3 occupied by 



two typical L* galaxies at z = 0, we obtain the UV LFs 
shown in Fig. [5] 

As seen from the same figure, at z ~ 6 and 2 ~ 7, both 
including/excluding the effects of dust on UV photons, the 
most luminous progenitors of the MW and M31 are bright 
enough to be amongst the faintest LBGs that have been 
detected to date. Furthermore, in the small range of overlap 
with the data at these redshifts, the slope and the amplitude 
of the theoretical UV LF is in reasonable agreement with 
the observations; this agreement is particularly striking at 
z ~ 7, where the deepest data are available, given that once 
the f c values have been calculated for any galaxy, there are 
no free parameters left to calculate the UV LF. However, at 
z ~ 8, even the predicted intrinsic UV LF of the MW/M31 
progenitors falls short of the current detection limits. 

Hence, our simulation suggests that, with the current 
detection limits, the most massive progenitors of local L* 
galaxies are luminous enough to only be visible at z ^ 7; 
above this redshift, the progenitors of local L* galaxies are 
too small to support enough star formation to be visible in 
the UV within the UDF12 data. Clearly our results thus 
also imply that high-redshift LBGs with Muv < — 18 are 
likely to be the progenitors of systems more massive than L* 
galaxies at 2 = 0. Finally, since the number of MW and M31 
progenitors increases dramatically with decreasing luminos- 
ity, it is clear that many such objects should be detectable 
with JWST, even at 2 ~ 8. 

3.3 UV Spectral Slopes 

We now explore the predicted colour excess, E(B — V) — 
— 2.51og(/ c )/11.082, produced by the dust attenuation in the 
simulated galaxies. We hav e used the SN extinction curve 
l|Bianchi fc Schneider! I2007T ) to infer the colour excess from 
the escape fraction of continuum photons for consistency 
with our assumption that SNII are the main dust produc- 
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Figure 3. The upper and lower panels show the colour excess, E(B — V), and the predicted UV spectral slope, 0, as a function of 
absolute UV magnitude for the progenitors of the MW (empty circles) and M31 (empty triangles) for z ~ 6 — 8 as marked above each 
panel. The vertical short (long) dashed lines show the detection limits for the HST (JWST). In panels (e) and (f), filled circles show 
the /3 slopes for individual galaxies, and filled sq uares show the average beta values in bins of one magnitude with the associated errors, 
inferred observationally bv lDunlop et alj l|2012bh . 



ers in these early galaxies. Further, the SN extinction curve 
has also been shown to success fully explain the obs erved 
properties of quasars at z ~ 6 (Maiolino ct al. 2006) and 
gamma-ray bursts at z = 6.3 (jStratta et al. 20071 ). 



Since we assume SNII to be the main dust producers, 
and the SNII rate is related to the total stellar mass that has 
been formed in stars, the more massive a galaxy is in M„, the 
larger its dust mass (and hence attenuation) is expected to 
be. This is exactly the trend seen from the upper 3 panels 
of Fig. [3] where the brightest progenitors show the largest 
colour excess such that E(B — V) ~ 0.02 while the faintest 
systems show negligibly-small values of E(B — V) ~ 0.001 at 
all the redshifts considered. The average value of E(B — V) 
increases from about 0.004 at z ~ 8 to 0.006 at z ~ 6, as 
progenitors increase in mass and build up their dust content. 

To obtain the observed spectral slopes, /?, for the pro- 
genitors of the MW and M31, we apply the SN extinction 
curve to the intrinsic spectrum of each progenitor, scaled 
according to the value of f c calculated for that progeni- 
tor. We use 30 values of /3 evenly spaced between —3 and 
to fit a line through the observed SED of each galaxy 
(as obtained using STARBURST99) over the wavelength range 
A rest = 1500 - 3000A sampled every 100 A; the value of 
/3 yielding the minimum \ 2 error is then chosen as the /3 
value for that particular galaxy. We find that /3 increases 
slightly with increasing mass at any given redshift as shown 
in the lower three panels of Fig. [3] for example, ft increases 



from about —2.64 to —2.35 as Muv increases from —13.0 to 
~ —17.5 at z ~ 6; this is due to increasing dust masses and 
attenuation values with increasing M* as discussed above. 

Including the inferred amounts of dust attenuation 
makes the observed /3 slightly redder than the intrin- 
sic values, such that the average value of ft is (/?) w 
(—2.54, —2.58, —2.6) for z m (6, 7, 8), as compared to an av- 
erage intrinsic ft value of (/?) — 2.64 at any of the redshifts 
considered. However, this relatively small change in average 
ft produced by dust is dominated by the very small amounts 
of dust attenuation in the large number of fainter progeni- 
tors. What is most striking in Fig. 3 is that the currently- 
detectable galaxies at these redshifts are predicted to have 
much more significant dust reddening, resulting in typical 
individual values f3 ~ —2.4. While this is still slightly bluer 
than the average value of = -2.02±0.29 (-1 .96 ±0.27) for 
Muv — —18 at z = 7(8) recently derived bv lDunlop et al.l 
|2012bl ), it is clear that the observable values predicted here 
lie well within the range of f3 = —4 to 1 (—3.5 to —0.8) ob- 
served for individual galaxi es at z ~ 7 (8) as show n in panels 
(e) and (f) of Fig. [3] (see also lDunlop et al.l2012al ). Moreover, 
since the volume sampled by the UDF12 at z ~ 6.5 — 7.5 is 
over an order of magnitude greater than that sampled here, 
the real galaxy samples at z ~ 7 and 8 inevitably contain 
many more massive galaxies than the most massive MW and 
M31 progenitors, and so the resulting sample-average values 
of f3 would be expected to be higher (i.e. redder) than the 
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UV slopes of even the most massive MW and M31 progeni- 
tors. 



4 SUMMARY AND CONCLUSIONS 

We have coupled state-of-the-art high-resolution gas- 
dynamical simulations of the local group run within the 
CLUES framework (that reproduce the mass, relative distance 
and velocity of the two MW and M31 haloes), with a dust 
model to explore the predicted observability and physical 
properties of the progenitors of the MW and M31 at z ~ 6, 7 
and 8 down to absolute UV magnitudes Muv < —13. Our 
main results can be summarized as follows. 

• At 2 ~ 6, 7 and 8, the simulated progenitors of the two 
L* galaxies at z — exhibit a three-dimensional relation 
between Af* — Z* — Muv, with the most massive progenitors 
being both the brightest in the UV, as well as the most metal 
rich. 

• The bright end of the predicted progenitor UV LFs at 
z ~ 6 (7) overlaps with the observed LFs for Muv — —17.25 
to —18.0 (—17.0 to —18.0). In this small range of overlap, 
both the slope and the amplitudes of the theoretical and 
observed LFs are in excellent agreement. However, at z ~ 8, 
even the intrinsic UV LF of the LG progenitors lies faint- 
ward of the observational limit at Muv — —17.25. 

• The above two results lead to a picture wherein at z ~ 
6 (7), analogues of the most massive progenitors of the MW 
(M31) are bright enough to be amongst the faintest LBGs 
that have been detected so far at these redshifts with Muv — 
-17.25 to -18.0 (-17.0 to -18.0). This implies that LBGs 
with Muv < —18 at z ~ 6 and 7 are likely the progenitors of 
systems more massive than L* galaxies at z — 0. At z ~ 8, 
even the most massive progenitors of the MW and M31 are 
intrinsically fainter than the current limit of Muv = —17.25 
implying that all LBGs observed so far at this redshift are 
the progenitors of super-L* galaxies at z = 0. 

• As a result of their relatively low stellar masses values 
at z ~ 6 (M* = 10 5,7-8 '°Mq), and even lower masses at 
higher redshift, the vast majority of the progenitors of the 
two LG galaxies have tiny colour excess values ranging from 
E(B - V) = 0.001 to 0.02. The average E(B - V) increases 
slightly from 0.004 at z ~ 8 to 0.006 at z ~ 6 as these 
galaxies build up their stellar and dust masses. 

• The dust attenuated UV spectral slopes become signifi- 
cantly redder with decreasing magnitude, from j3 ~ —2.64 to 
—2.35 for Muv — —13 to —17.5 at z ~ 6 as a result of the 
larger dust masses in the most massive (luminous) galax- 
ies. The brightest progenitors at z ~ 7 with Muv = —17.4 
and —17.8 exhibit /3 values of —2.5 and —2.4 respectively. 
These values are slightly bluer than the average value of 
/9 = -2.02 ± 0.2 9 inferred for M uv = -17.5 LBGs at the 
same redshift bv lDunlop et al.l (|2012bh . but lie well within 
the observed scatter that ranges between f3 = —4 to 1. This 
is broadly as expected, given that the most massive progen- 
itors of the MW and M31 at these redshifts are comparable 
to the least luminous gal axies in the UDF12 samples studied 
bv lDunlop et al.1 (|2012tJ) . 
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